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The aim of the present paper is the study of the entropic elasticity of the dsDNA molecule, 
having a cristallographic length L of the order of 10 to 30 persistence lengths A, when it is subject 
to spatial obstructions. We have not tried to obtain the single molecule partition function by solving 
a Schodringer-like equation. We prefer to stay within a discretized version of the WLC model with 
an added one-monomer potential, simulating the spatial constraints. We derived directly from 
the discretized Boltzmann formula the transfer matrix connecting the partition functions relative 
to adjacent "effective monomers". We have plugged adequate Dirac (^-functions in the functional 
integral to ensure that the monomer coordinate and the tangent vector are independent variables. 
The partition function is, then, given by an iterative process which is both numerically efficient and 
physically transparent. As a test of our discretized approach, we have studied two configurations 
involving a dsDNA molecule confined between a pair of parallel plates. One molecule end is anchored 
to one plate by a biochemical bond. A stretching force F, normal to the plates, is pulling away the 
other end. In the first case, the cristallographic length L is smaller than the two-plate distance Lq. 
The molecule feels, then, only the anchoring barrier effect. The predicted elongation-versus-force 
curve, is pushed upward with respect to the WLC model result. This effect is the most spectacular in 
the low force regime. For large forces say, Fh% g h = 5 ks T/A, the elongation versus L is very well fitted 
by a straight line with a slope given by the standard WLC model and a constant term ~ 1.2^4. In 
the second case, L takes values up to L max = 1.5 Lq. With a stretching force still equal to Fhigh, the 
standard WLC model predicts that the molecule cannot fit within the plates when L > L, — 1.29Lo- 
We have studied the evolution of the elongation derivative with respect to L, together with the mean 
square free-end fluctuations along the force. They both exibit a sharp decrease when L > Lq. We 
present a semi-qantitative argument suggesting that the terminal segment involving 20 % of the 
internal monomers flattens against the repulsive barrier when L — > L max . In conclusion, we suggest 
extensions of the present work, relevant to the analysis of micromanipulation experiments. Finally, 
we have gathered into the Appendix formal developments, leading to a precise relation between the 
transfer matrix and the Hamiltonian methods for the study of spatially constrained dsDNA. 

PACS numbers: 87.15.By, 61.41+e 

Introduction 

In the last ten years Single Molecule Biophysics has become a very active field of research. Among the recently 
explored areas, one finds the observation, at the one-molecule level, of the biochemical interactions of the double 
strand DNA (dsDNA) with the various proteins involved in the duplication process (for two recent reviews see, for 
instance, the references 0,0)- The protein-DNA interaction is detected by the observation, in real time, of the 
variations of the dsDNA elongation under the action of a fixed stretchingforce. It was clearly of interest to have a 
good physical understanding of the dsDNA elasticity measurements 13, 0,13 ■ I n absence of DNA supercoiling, this is 
provided by the so called "Worm Like Chain" (WLC) model 0, S Hi which gives a good description of the dsDNA 
entropic bending elasticity within a wide range of force, from few hundredths to few tens of picoNewton. 

In practice, the persistence length of a double helix is about five times the typical length resolution. It is, then, 
legitimate to use a rectifiable curve to represent the coarse-grained dsDNA chain. The basic hypothesis of the WLC 
model is to assume that the elastic-energy linear density is inversely proportional to the square of the chain-curvature 
radius. In the usual formulation of the WLC model, the relevant Statistical Mechanics variable is the tangent vector 
t(s) = dr/ds, where r(s) is the effective monomer coordinate and ds the rectifiable chain length element. The clastic 
molecular-chain energy, expressed in thermal unit ks T, is given by the following line integral: 

E WLC =J o L d S ^(^ir-i-t( S) y (i) 

where L is the crystallographic length of the polymer. A is the persistence length and (dt(s) / ds) 2 is the inverse 
square of the curvature radius of the coarse-grained chain. F = ks Tf stands for the stretching force, applied to the 
extremity of the chain. It follows from its very definition that t is a unit vector. The partition function Z is given by 
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the functional Boltzmann integral: 

Z= fv [0]D [<fi}exp(-E WLC ) , (2) 

where the integration has to be performed over all the paths joining two points on the unit sphere. Exploiting the 
analogy with Feynman formulation of Quantum Mechanics (QM), the result of the functional integration is given in 
terms of the Hamiltonian operator [l8| : 

iW = -^V t 2 -f t, (3) 

where V 2 is the Laplacian on the unit sphere. In the limit where the L > 4, the elongation of the molecule under 
the action of a stretching force F = ksT f directed along the z — axis is given by: 

(z(L)) = -L°ffi, (4) 

where Eo(f), the ground state eigen-energy of the Hamiltonian Hwlc, can be obtained with a very good precision 
by solving numerically an ordinary differential equation with appropriate boundary conditions Q . 

A. The problem of the Spatial Constraints in the WLC Model. 

However, this approach becomes ackward if one wishes to impose physical constraints involving the monomer space 
coordinate. For instance, one may wish to confine the molecule within a certain region of space. This constraint can 
be formulated in terms of a potential, V(r(s)), acting on each monomer of coordinate r(s). If the unitary tangent 
vector t(s) is the sole dynamical variable, one has to write r(s) = J Q S t(s')ds' and the potential energy to be added to 
Ewlc has the non-local form: 

AE^ LC = £dsv(J\(s')dA . (5) 

This situation becomes worse if one studies self-avoiding effects described by a monomer- monomer repulsive potential, 
y(r( Sl ),r( S2 )): 



AEwlc =[ ds i [ ds2 V ([ l t(ai)d«i , / 2 t(s' 2 )ds' 2 
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(6) 



The obvious thing to avoid this problem for the monomer potential energy is to use the monomer coordinate r(s) 
as a dynamical variable. It is convenient to consider a whole family of polymer models defined by the elastic-energy 
linear density: 

£( S )=£o(T 2 ) + ±Ar 2 -f-r + V(r), (7) 

where the WLC model can be recovered by an appropriate choice of £o(r 2 ). But there is, clearly, a price to be paid : 
second-order derivatives appear now in the elastic-energy linear density, via the curvature term ^ A f 2 . The mapping 
onto an Euclidian Quantum Mechanics problem is no longer as evident as it was before. 

B. A very brief review of previous studies on confined semiflexible polymers. 

Numerous authors have addressed themselves to the problem of fin ding an Hamiltonian describing semi-flexible 
polymers subjected to spatial constraints [T^. ll3L HI . IT5I TT6I IrH fl8l IT^ . l20t I2TL |22 | . The basic idea was to formulate 
the problem in such a way that both the coordinate r and the "velocity" v = r appear as independent dynamical 
variables. Various mathematical techniques have been used to arrive to the Hamiltonian : 

H = -^V v 2 + £ (v 2 ) + v • (V r - f) + V(r) . (8) 
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The case of the WLC model is obtained by choosing £o(v 2 ) *- v 2 g ^ and going to the limit Sb — > 0. It leads to the 
spatially constrained WLC Hamiltonian: 

Hscwlc = H WLC + t-V r + V(r). (9) 

A group of authors [T3. IT^L IbH . ITHl ITfiL IT^ | have focused their analysis on the case of semiflexible polymers confined 
within a cylindrical tube having a radius much smaller than the persistence length. It corresponds to the choice 
£o(v 2 ) = 1 with the running variable s being identified with the monomer-coordinate component along the axis of 
the tube. This can be viewed as the statistical physics analog of the Monge approximation for elastic rods. This 
model has been applied recently to interesting physical problems: the first one is the theoretical study of the flow of 
semiflexible polymers across cylindrical pores [T|j, a second one is the analysis of the unbinding transition between 
semiflexible polymers and directed polymers acting as one-dimensional attractive systems. |17| . 

The fully three-dimensional WLC Hamiltonian, involving both the coordinate r and the "velocity" v — t, appears 
in references fl8l Il9ll20ll2lll2^ | . One of the most significant applications, which bears some resemblance with the work 
presented here, is found in reference p3 |. where the authors analyse a symmetric interface between two immiscible, 
semiflexible polymers. 



C. Organization and Synopsis of the Paper 

The aim of the present paper is the study of the spatial obstructions which are present in most micromanipulation 
experiments, involving a single dsDNA, immersed in a liquid thermal bath and possibly interacting with a single 
protein. These obstructive effects are expected to be significant in recent el ong ation experiments, involving molecular 
segments with a cristallographic length L of about ten persistence lengths |27| . 

We have not tried to obtain the partition function of a ds-DNA molecule, subject to spatial constraints, by solving 
the partial differential equation associated with a Schrodinger-like problem, written with an imaginary time equal 
to —i s. In Section I, we rather stay within a discretized version of the WLC Model involving N links (or effective 
monomers) of length b = L /A. The partition function Zn is obtained by an iteration procedure involving the "transfer 
matrix" connecting adjacent monomers j24[. The "velocities" v„, relative to the N effective monomers, are introduced 
in the discretized Boltzmann formula via appropriate Dirac ^-functions. The WLC model is obtained by choosing: 
£o( v n) = b/(2Sb 2 ) (v n — l) 2 and by performing the integrals over the moduli v n = |v„| in the limit Sb/b <C 1, so that 
the tips of v ra are restricted to the unit sphere. After few manipulations, we arrive to a recurrence relation for the 
partition functions of molecules having their lengths L which increase by unit of b: 

Z n+ i(r„ + i,t„ + i) =exp(-&V(r n+ i)) J d 2 n(t n )T WLC (t n +i | t„) Z n (r n+1 -6t n+ i,t n ). (10) 

Here Twlc{^u+i |t n ) is the transfer matrix relevant for the unconstrained WLC model. The above iteration pro- 
cedure, which will be our basic tool for the study of the entropic elasticity of dsDNA subject to spatial constraints, 
exhibits a suggestive connection with a Markovian random walk in three dimensions. 

The Section II is devoted to semirealistic applications of the above formalism. Our testing ground is the study of the 
entropic elasticity of a single dsDNA molecule confined between two parallel plates. One molecular end is anchored 
to one plate by a biochemical binding. A stretching force, normal to the plates, is pulling away the free molecular 
end from the anchoring plate. As long as we are mainly interested in the elongation of the ds-DNA molecule, we can 
exploit the invariance under translations parallel to the plates and rotations around the stretching force to write a 
recurrence formula involving only longitudinal variables: 

Z n+1 (z n+1 ,8 n+1 ) = exp (~bV(z n+1 )) / d(cos9 n )T WLC (9 n +i, 9 n , f)Z n (z n+1 - b cos0 n+ i, 9 n ), 

Jo 

where z and cos 9 are respectively the components along the force direction of the coordinate r and the unitary velocity 
vector t. 

For the sake of simplicity, we ignore, in this preliminary analysis, possible spatial obstructions associated with the 
stretching devices (magnetic and optical tweezers), but, in Section III, we suggest a practical way to take them into 
account. Two configurations are studied. 

In the first one, the cristallographic length L is supposed to be shorter than the two-plate distance Lq, so that 
the entropic elasticity is affected only by the anchoring plate barrier. For a relatively short molecule (L/A=12), we 
compute the elongation-versus-force curve, which is expected to be pushed upward with respect to the unconstrained 
WLC model predictions, notably for the zero-force case. In the case of relatively high force F, i.e a — F A/(ks T) = 5, 
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the elongation ( z(L) ) versus L is very well fitted, within the interval 2 A < L < 10 A, by a straight line with a slope = 
0.771 and a constant term = 1.22 A. The slope turns out to be very close to the value predicted by the unconstrained 
WLC model in the limit A <C L, namely, 0.775. If the second plate is pushed further away, the predictions of the two 
models for ( z(L) ) /L will coincide when A <C L up to a 0.5% correction, which could be due to our use of a discretized 
version of the WLC model. 

In the second configuration, L is allowed to take values up to 1.5 Lq, with the same stretching force value as above: 
a = F Aj [ks T) = 5. A molecule, elongated according to the unconstrained WLC model cannot fit within the 
plates when its cristallographic length L is larger than the critical value = 1.29Lo, associated with the elongation 
( z(L*) ) = Lq. We have studied the evolution of the terminal monomer statistics when L varies within the interval: 
0.6 Lq < L < 1.5 Lq, using as representative quantities the elongation derivative d ( z{L) ) / dL and the mean square of 
the free-end fluctuations along the stretching force direction, i.e. Az 2 (L) = ((z(L) — (z(L))) 2 ). When L < L these 
two quantities follow rather closely the predictions of the unconstrained WLC model but they both start a rather 
sharp decrease when L > Lq. The elongation (z(L)) is no longer an extensive physical quantity an goes slowly to 
Lq. As to the fluctuations Az 2 (L), they are reduced by a factor 10 with respect to the unconstrained WLC model 
prediction. 

In the case L = L max = 1.5 Lq, we give a semi-quantitative analysis of the internal monomers statistics, when the 
monomer number varies within the interval L*/b < n < L max /b. The characteristic feature of our model lies in the 
fact that all the monomers are confined between the two plates. We call this type of spatial constraint an internal 
confinement (IC) in contrast with the external confinement (EC) describing a situation where the confining V{z) is 
acting only upon the terminal monomer. A rather good physical understanding of what is going on near the repulsive 
barrier has been obtained within a Gaussian model, supplemented by a square- well potential with a depth 3> kgT. 
Choosing a configuration with the same value of L„ as before, we find, without too much surprise, that in the IC 
model the terminal segment involving 20 % of the internal monomers flattens against the repulsive barrier, while in 
the EC model all the monomers, except the terminal one, can move rather freely across the repulsive barrier. 

In Section III, we suggest possible applications or extensions of the work presented in the present paper; let us 
mention the two of them which concern directly the micro-manipulation cxpriments: 

• We propose a simple procedure involving a two-plate confinement model which may lead to an estimate of the 
spatial obstruction effects associated with the magnetic tweezer. 

• It is relatively straightforward to generalize the spatially constrained WLC model, within its transfer matrix 
formulation, to the RLC model 0, El El, which incorporates both bending and twisting rigidities. The 
anchoring barrier effects are expected to be significant when the reduced supercoiling parameter a is above the 
threshold - at fixed force - for the creation of plectonem configurations. The stretching potential energy vanishes 
for such structures, allowing them to wander irrespective of the sign of f ■ t. 

The Appendix is devoted to a comparison of two possible approaches to the elastic entropy elasticity of dsDNA 
subject to spatial constraints. The first one, used in the present paper, is the transfer matrix method which leads to 
an iterative construction of the par tition funct ion, within a discretized version of the constrained WLC model. The 
second approach [1J, [lj, [lj, [l5|, Ha, Hll 0, Ha |2(J, EDI is based upon the solution of a Schrodinger-like equation, 
written with an imaginary time variable —i s involving the Hamiltonian given by eq. J^J. To achieve our purpose, we 
have found convenient to use the auxiliary variable method in order to eliminate the second order derivative in the 
clastic energy density given by eq. (|12fl . It is, then, a rather straightforward affair to obtain the transfer operator 
T- associated with the transfer matrix-, from which one derives the Hamiltonian H by taking the continuous limit 
b — > 0. Making use of the leeway inherent to any discretization procedure, one can obtain a symmetric version of T, 
which coincides with the exact evolution operator exp—bH, up to corrections of the order of b . 

It is not too difficult to elucidate the physical interpretation of the auxiliary variable u: we have proved that the 
conjugate momentum p u appearing in H is just the "velocity" v. Performing an appropriate change of basis, one 
obtains immediately an expression of H, identical to the r.h.s. of eq.(|HJ). With the proper choice of £q{y 2 ), one then 
arrives to the spatially constrained WLC Hamiltonian, Hscwlc, given by eq.© [llOQIIlElllHIlilil 
I2ll l22j| . Finally, as an internal check, we derive a symmetric transfer matrix from the transfer operator T, associated 
with the Hamiltonian Hscwlc- By making suitable approximations, we recover the transfer matrix, physically more 
transparent, which is derived directly from the Boltzmann formula in Section I. 



5 



I. A TRANSFER MATRIX APPROACH TO STRETCHED DS-DNA SUBJECT TO SPATIAL 

CONSTRAINTS 

We wish, first, to consider the statistical properties of a class of polymer models described by the linear elastic- 
energy density £(s) depending upon the monomer coordinate r(s), together with its first and second order derivatives, 

r = 4| and r = 4-z (r) . The variable s with < s < L results from a coarse graining of the molecular chain, having a 
cristallographic length L. Note that s does not coincide necessarily with the arc-length s of a rectifiable curve. 

The partition function for fixed free-ends polymers is given by a functional integral involving the standard Boltzmann 
stastitical weights : 

Z = yp[r]exp^-^ £(s)ds^j , (11) 

£(s)=£ (i 2 ) + ^Ar 2 -{v + V(v), (12) 

where the functional integration goes over all the paths joining the free ends of the chain of fixed coordinates r(0) 
and r(L). 

We stress that, here, the potential V(r), simulating the spatial constraints, is clearly acting upon all the monomers 
of the molecular chain, contrary to confining models where the potential V(r) is acting, only, upon the terminal 
monomer, e.g. via an attached bead, In this latter case, constraints can be trivially implemented by adding in the 
elastic density <|12|) the "velocity" dependent contribution r ■ W r V(r). 

To compute the partition function we have to resort to a discretization of the variable s: s n = nb. The molecular 
chain is then represented by N elementary links or effective monomers with N = L/b. Assuming that the effective 
monomer length b is much smaller than the persistence length A, we can write: 

I'u i",,. i .. r„ — 2r n _i + r„_2 



b " b 2 

It is then convenient to introduce the discretized form of the elastic-energy density £disc(n) obtained by plugging the 
formulas Ijl3|l into the right hand side (r.h.s.) of equation l|12|l: 



w») = £o ( r ;-' i ) 2 j + \ a p ""-^~-^ )a _ f . r + v{Yn) _ (14) 

The partition function is then given as an integral of a product of N — 1 monomer-coordinate-r„ functions: 

„ N-l N 

Z= / [] d z r n {\[e^{-b£ dlsc {n))}Z m {^o), (15) 

J n=l n=2 

where the monomers coordinates ro and are kept fixed. The starting partition function is taken to be : 

Z iB (n,r ) = exp (-b£ ffi^°)A - bV(rx) \ . (16) 

From the above expression of the partition function one could write down easily a transfer matrix connecting the 
probability distributions relative to pairs of adjacent monomers: 

Z n+1 (r n+1 ,r n ) = J d 3 r„_id 3 r„_ 2 exp(-&£ disc (n + 1) - b£ disc (n)) Z„_i(r„_i, r„_ 2 ). (17) 

This approach turns out to be rather awkward in the particular case of the WLC model but may be useful for 
other models. The standard trick is to introduce the "velocities " v„ = r„ = (r„ — r„_i)/& as new variables in the 
discretized functional integral l|15[) via the trivial identity: 



1 = J d 3 v n * 3 (v n - r " b r - 1 ). 



Making some obvious manipulations, one gets a new expression for the partition function as an integral over paths 
joining two elements in the 6-dimension space obtained by taking the direct sum of the coordinate r and velocity v 
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spaces: 

, N-l N 



Z N (r N ,v N ) = f l[d s v n d 3 v n {f[ 5 3 (v n - rn r "~ 1 ) 

« 1 n 



n=2 

cxp -6 ^o(v^) + ( V "Y"" X ) 2 - f • v„ + V(T n )\ } Zi(n, Vl ). (18) 

For practical reasons, the coordinate ro will be assumed from now on to be a random variable with the probability 
distribution -P(r ). Then, Zi(ri,Vi) is given by: 

Zi(n, vi) = exp-6 (f (v?) - f • vi + F(n)) P (n - 6vi). (19) 

If one performs explicitly the integration upon the v n variables in the r.h.s of equation JTSJ, one readily recovers the 
r.h.s of equation (|15fl with Z in {vi, r ) replaced by its average over r . Only nearest-neighbour monomers with given 
coordinate and velocity are now connected. It is then possible to write down a recurrence relation between adjacent 
intermediate partition functions Z n (r n ,v n ) relative to chains having a cristallographic length s n = bn: 

Z„ + l(r n+1 ,v n+ i) = J J d 3 r n d 3 v n T(r n+1 ,v n+1 | r„, v„) Z„(r„, v„). (20) 

The transfer matrix T(r n+ i, v n +i | r„, u„) is easily read off from the r.h.s of equation l|18|) : 

T(r„ + i, v„ + i | r„, v„) = 5 3 (v n+ i - r " +1 r " ) exp -b ( V(r n+1 )) x 

exp -b (W 2 +1 ) + i A ( V " +1 6 ' V " ) 2 - f • v n+1 ) . (21) 

Performing the integration upon r„, we arrive to the the final form of the recurrence relation: 



Z n+ i(r n+1 ,v n+1 ) = exp-b (£ (v„+i) 2 + V(r n+1 )) Jd 



3 v« x 



cxp- A(v n+ i - v„) 2 - bf ■ v„ + i^ Z n (r n+X - bv n+ i ,v n ). (22) 

In order to apply the above formula to the WLC model with spatial constraints, we have to choose an adequate form 
for the function £o(v 2 ). The solution is rather simple: one introduces the small length 8b such that 5b/b <C 1 and 
one takes for £o(v 2 ) the following expression: 

fo(v - )=6 H^- (23) 

Indeed, with this choice, we write the volume element in the velocity space as: J d 3 v n = Vndv n d 2 n(t n ), where 
d 2 Q(t n ) is the infinitesimal solid angle around the unit vector t„ taken along v„, written as v n t n . Then, in the limit 
Sb/b — > , exp — 6 (£o(v 2 )) reduces, up to a numerical constant, to the 5-function: S(v n — 1). As a consequence, 
integrating over v n within the same limit, leads to an expression of the r.h.s of equation l|18(l readily obtained by 
making the replacement: 

d 3 v„ — > d 2 fl(t n ) , v„ — > t n . 

The transfer matrix describing stretched dsDNA within the WLC model with spatial constraints is then simply 
obtained from ea. (|22|l using the same replacement rules: 



Z n+1 (r n+1 ,t n+1 ) = exp -b (V(r n+ i)) J d 2 tt(t 7 



cxp- ^(t„+i-t„) 2 -f • (t n+ i + t„)/2j Z n (r n+ i -6t n+ i,t„). (24) 

To make easier the comparison with the standard WLC model, we have replaced in the above equation f • t Jl+1 
by f • (t n +i + t„)/2. Such a modification amounts to the replacement of a rectagular discretized integration by a 
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trapezoidal one. What appears now in the r.h.s of the recurrence relation is the transfer matrix relative to the standard 
WLC model: 

T WLC (t n+1 1 1„) = exp (- A(t " + 2 1 & ~ t " )2 + \ f ■ (Wi + t„) ) . (25) 

In the Appendix, we are going to show, within a more general context, that T\YLc(^n+i | t n ) coincides, up to corrections 
of the order of (b/A) 3 , with the exact transfer operator exp (—bHwLcj with Hy/LC = ~TA^t — f • t. 

We are, now, ready to give the partition-function recurrence relation which will be our basic tool for the study of 
the entropic elasticity of dsDNA subject to spatial constraints: 

Z n+1 (r n+1 ,t n+1 ) = exp (-6V(r n+ i)) J d 2 n(t n )T WLC {t n+ i \ t„) Z n (r n+1 - bt n+1 , t„) . (26) 

The iterative construction of the partition function Z n (r n ,t n ) has a transparent physical meaning. Indeed, it has a 
suggestive interpretation in terms of a Markovian random walk model in three dimensions: the (n + l) th step is given 
by bt n+ i; its length is b and its direction, defined by the unitary vector t n+ i, is correlated to that of the previous step 
t„ through the correlation function C(t„ + i, t„) oc Twlc{^u+i | tn)- Within our unit conventions, the "potential" V(r) 
stands for a potential energy by unit molecular length, written in thermal units, so that AW n +i = kBTbV{v n+ \) 
represents the amount of energy to be exchanged with the thermal bath by performing the (n + l) th step. The 
exponential in front of the r.h.s of equation l|2fci[l is just the associated Boltzmann factor. 

II. THE STRETCHING OF A DS-DNA MOLECULE CONFINED BETWEEN TWO PARALLEL 

PLATES. 

In this section, we apply our iterative version ea. ^26|) of the spatially constrained WLC model to the confinement of 
a single dsDNA molecule between two parallel plates - called hereafter Plate 1 and Plate 2 - separated by a distance 
Lq. One molecular end is anchored upon Plate 1. The other end is attached, for instance, to a magnetic bead pulled by 
a magnetic tweezer. We do not pretend, here, to describe a fully realistic situation since, for the sake of simplicity, we 
ignore the spatial constraints associated with the finite dimensions of the magnetic bead. In Section III of this paper, 
we shall suggest a self-consistent procedure, within the two-plate model, which will lead to an estimate of the bead 
obstruction effect. The stretching force is assumed to be normal to the two plates and is pulling away the molecular 
end from the anchoring plate. The dimensions of the plates normal to the force, measured from the anchoring point, 
are assumed to be much larger than Lq and will be considered as infinite. Two physical configurations have been 
considered: 

1. Lq is supposed to be larger than the cristallographic length L. (In our explicit computations we have taken 
Lq — 40 A and L = 12 A.) It means that even for large stretching forces the molecular end stays far away from 
Plate 2. This will allow us to study the deviations of the elongation from the WLC model prediction due to the 
presence of the anchoring plate (Plate 1), provided we ignore the spatial obstruction of the pulling device. 

2. The cristallographic length L is larger than Lq so that, for large forces, the fully stretched molecule does not 
fit inside the plates. We have studied, within the WLC model, how the molecule adapts itself to a confined 
situation where the elongation predicted by the unconstrained WLC Model reaches values significantly larger 
than Lq. 

A. Basic formulae for the two-plate confinement configurations. 

Our confinement configurations, defined by the two confining plates and the stretching force, are invariant, first, 
under rotations around the z-axis defined by the anchoring point and the stretching force direction, second, under the 
translations within the x, y plane. We shall consider a molecular ensemble where the anchoring monomer is uniformely 
distributed over the x, y plane. The iteration procedure can then be organized in such a way that at each step the 
partition function Z n (r n ,t„) depends only upon the two independent variables: z n and cos# n , which are respectively 
the components upon the z-axis of the monomer coordinate r„ and the tangent vector t„. If this condition is satisfied 
at the step n, Z n+ i(z n+ i, 9 n +i) is then easily obtained from an azymuthal average around the z-axis of the two sides 
of eq.jgg): 
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FIG. 1: The two-plate confining potential V(z) (dotted line) with the associated Boltzmann factor exp — bV(z) (full line). We 
have taken as elementary link length b — A/10 . The Boltzmann factor is given explicitly by the product of two rounded-off 
step functions Q(z, Az) Q(Lq — z, Az) with a smoothing length Az — A and a two-plate distance Lo = 4CL4. 



Z n +i(z n+ i,0 n+ i) = exp (-6 V(z n+ i)) / d(cos9 n )TwLc(0n+i,0 n , f) x 

Jo 

Z n (z n +i-b cos 6 n+ x,6 n ) (27) 



The azymuthal-averaged WLC transfer matrix can be obtained by taking the zero torque limit of the corresponding 
expression given in ref.[ll| for the supercoiled DNA case: 

A b f 

Twlc{Qi, 02, /) = exp-{— (1 - cos0i cos0a) + — (cos#i + cosfe)} 

A sin 6i sin 6 2 , , ocA 
x Io( ~ b ) • (28) 

In order to build a potential V(z) simulating the two confining plates and satisfying a smoothness condition at the 
persistence-length scale, it is convenient to introduce the rounded-off step function: 

6(z,Az) = i + ierf(z/Az), (29) 

where erf(cc) is the "error" function : cxp(— t 2 ) dt and Az the smoothing length assumed to be ~ A. Instead 

of the potential V{z), it is more convenient to write directly the Boltzman factor : 

Boltz(z, L , Az) = exp (-b V{z)) = Q(z, Az) 0(L Q - z, Az) . (30) 

Its variation with respect to the reduced variable z/A is represented by the solid curve of Figure 1 for Az = A. The 
dotted curve gives the variation of the external potential associated with an elementary link bV(z n ). 

The next point to be specified is the inital condition of our recurrence procedure. A detailed analysis of the 
anchoring mechanism at the microscopic level being beyond the scope of this paper, we have adopted a prescription 
which incorporates in a simple way some features of the actual physical situation. We choose as the molecular chain 
origin (n = 0) the end of the initial strand of length ~ A sticking out from the anchoring plate. The initial partition 
function Zq(zq,9o) is written as follows : 

Z (zo, ) = - t- exp (-h Za /a) 2 ) . (31) 
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FIG. 2: Illustration of the effects of the anchoring plate barrier upon the elongation-versus-force curve. The upper curve 
(full line) gives the prediction of the spatially-constained WLC model when L/A = 12. The other plate is lying at a distance 
Lo = 10 L/3 and, as a consequence, has no effect upon the elongation. The fact that the elongation curve is pushed upward with 
respect to the unconstrained WLC prediction(dotted curve) has a simple qualitative explanation: because of the anchoring 
plate barrier the z < half space is not accessible to the nucleotides, whatever the stretching force. 



A look at the potential V(z) on Fig. 1 shows that zo can take negative as well as positive values within the range 
(—A, A), so that a centered Gaussian curve looks quite realistic if a ~ A. Similarly the tangent-vector z-axis projection, 
cos#0j hes within the interval (—1, 1), so that a cos#o uniform distribution seems to be an acceptable guess. We note, 
finally, that our choice of Zq(zq,9q) satisfies the smoothness condition assumed in the previous section. The n = 1 
partition function Z\(z\,9\) is then readily obtained by plugging Zq(zq,9q) in the r.h.s. of the recurrence relation 

ED 

We would like to give, now, a few indications about the numerical methods we have used to run the transfer matrix 
iteration process. 

In order to perform the integral over 9 ni we use the following discretization procedure. We divide the variation 
interval < 9 n < n into n s segments ^~ 1 ^ 7r < 9 n < — . The integral over each segment is done with the stan- 
dard Gauss method involving n g abscissae and n g attached weights. The integral over the full 9 n interval is then 
approximately given by a discrete weighted sum over d = n s n g points: 

/ f(9 n ) sm6 n d0„ =y^Wj sm0if(6i) . 

J ° i=l 

The transfer-matrix iteration involves a z variable translation z n — > z n+ \ — bcos9 n+ i to be performed upon 
Z n (z n ,9 n ) for each value 9 n = 9i. This is achieved by building at each step the interpolating function Z n ^ n t{z,6i) 
associated with the array: {zi = z m i n + I b, Z n (zi, 9i)} with 1 < I < n max . An appropriate choice of n max leads to a 
physically relevant sample of the monomer coordinates throughout the iteration process. 

For more extensive computations than those presented in this paper, one should consider a potentially more ef- 
ficient method involving the Fast Fourier Transform (FFT) algorithm. It is based upon the remark that a vari- 
able translation performed upon a given function reduces to a phase shift upon its Fourier transform: Z{p 1 9i) — > 
exp (— ibcos9i) Z(p, 9i). An inverse FFT will then be required to perform the multiplication by the Boltzmann factor. 



B. Modification of the elongation-versus-force curve induced by the anchoring plate barrier. 

We would like to present, here, the results of a numerical simulation based upon the iteration process built up from 
the recurrence relation given in ea. l|27|) . Our aim was to study the barrier effect of the anchoring plate (Plate 1.) upon 
the elongation-versus-force curve. This corresponds to the configuration 1, introduced previously, where the distance 
between the plates Lq is larger than the cristallographic length L. The values used in our simulation are Lq = 40 A 
and L = 12 A. It is, then, clear that Plate-2 plays no role since we have decided to ignore the eventual presence of a 
magnetic bead attached to the molecule free end. At each step n we store the normalized z„-probability distribution, 
P n (z n , a), where a = F Aj{k,B T) is the reduced force parameter. We have plotted in Figure 2 the relative molecular 
elongation ( z(L) )/L = J dzjy z^ Pn(zn, Oi)/L versus a, where N = L/b. The dotted curve represents the prediction 
of the standard WLC model which assumes that the anchoring device reduces to a point. 

The full curve corresponds to the elongation-versus-force obtained in our simulation within a spatially constrained 
WLC model; it exhibits clearly the anchoring-plate barrier effect which leads to a relative elongation of 34% in the zero 
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FIG. 3: Linear fit to a set of ( z(L) ) values, coming out from our iteration procedure in the case a = 5, when L increases from 
2 A to 10 A. The slope given by the fit, 0.771 is in good agreement with the prediction of the standard WLC model, 0775, valid 
in the limit A/L < 1. 



force limit. For higher forces, the two curves become approximately parallel with an offset of about 10%. It is possible 
to get from our simulation an expansion in powers of -J A/L at fixed force, giving the upward elongation displacement 
Ab ar (z/ L) induced by the anchoring plate barrier. As an illustration, we give the result for two typical forces 

corresponding to a — and a = 1. We found that the second-order expansion, Af, ar (z/L) = ai(a) \J^+ 0.2(a) j;, gives 

a very good fit to our simulation data within the range: 5A<L< 30 A. For a = we get: ai(0) = 1.00 , a 2 (0) = 0.65. 
The large value of ai(0) can be interpreted in terms of an hemispheric molecular cluster around the anchoring point 
with a radius growing like y/N. The results look rather different when a = 1: ai(l) = —.06 , 02(1) = 1.57. The strong 

decrease of the coefficient - confirmed by the a = 2 results - seems to suggest that the anchoring plate barrier 

affects only a dsDNA segment having a fixed length ~ A when a > 1. 

When one looks at the curves of Fig. 2, one may wonder if they are going to meet in the limit of high force, as 
they should, since (z(L))/L goes to one for the two cases. Before trying to push the present computation to larger 
values of a, one has to keep in mind two things, one, the dotted curve corresponds to the limit while we 

are working with A/L ~ 0.1, second, the dotted curve goes rather slowly to the limit a — > 00: indeed, (z(L))/L 
=0.95 for a ~ 100. Our computation will simply not work under these two extreme conditions. There is fortunately 
another way to make a meaningful comparison of the two models: it is to look at the derivatives d<K z ^ - instead 
of (z(L))/L. In the unconstrained WLC model computations these two quantities arc identical but this is not in 
general true in presence of spatial constraints. We have plotted in Fig. 3 the values ( z(L) ) obtained by our iteration 
procedure in the case a = 5. The results correspond to the big dots and cover the range of L values: 2 A < L < 10 A. 
We have performed a linear fit to the data, (z(L))fu — 0.771L + 1.22A, which, as it is apparent on Fig. 3, works 
beautifully. The fitted slope, 0.771, is very close to that given by the unconstrained WLC model, namely, 0.775. We 
also note that the constant term 1.22 A has the same order of magnitude as the one obtained previously for a = 1. 
This strongly suggests that, if the Plate-2 is pushed to infinity, the two models will give identical predictions, for the 
elongation (z(L))/L in the limit A<L. We have ignored the 0.5% difference, which may be attributed to our use 
of a discretized model wih b/A = 0.1. 



C. Statistical behaviour of stretched molecules having cristallographic lengths larger than the two-plate 

distance. 

This section is devoted to the analysis of the results of the transfer-matrix iteration for molecular chains having 
cristallographic lengths within the interval 20A < L < 30A when they are confined between two plates separated by 
a distance Lq — 20A The stretching force is specified by taking the value a = F A/(ks T) = 5 for the reduced force 
parameter. The unconstrained WLC model would lead to a relative elongation (z(L)) /L — 0.775, to be compared with 
the ratio L /L — 2/3 ~ 0.667. As a consequence, molecules having a cristallographic length L > L* = i /0.775 = 
1.29 Lq = 25.8 A, could not fit within the two plates if stretched according to the WLC model prediction. 
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FIG. 4: 3-D plot illustrating the probability distribution of the coordinate z of the molecular end, P(z,L), when the cristallo- 
graphic length L varies from 0.9 Lo to 1.5 Lo with Lo = 20 A. In our iterative computation, we have taken a stretching force 
such that the elongation predicted by the unconstrained WLC model is { z(L) ) = I.I6L0 when L takes its maximum value. It 
is then clear that such a molecule, stretched according to the usual WLC model, could not fit between the plates. When L 
approaches its maximum value 1.5 Lo, the elongation ( z(L)) > is no longer increasing like L but tends to its maximum possible 
value Lo; at the same time the slope of the hill exhibited by P(z,L) is becoming steeper. This indicates a decrease of the 
longitudinal fluctuations of the terminal monomer. 



Our iteration procedure provides us with the the terminal-monomer distributions P{z, L) for various molecular 
lengths L < 30 A. At each step we store the probability P n (z(n)) and by identifying L = nb we build the discrete 
set: P(z,nb) = P n (z). The probability distribution P(z,L) is then obtained by interpolation. We have displayed in 
Fig. 4 a 3-D plot of P(z,L). It allows us to follow in a continuous way how the DNA molecule manages to satisfy 
the two-plate space constraints. When 18 < L < 20 the molecule does not feel yet the Plate-2 barrier and the hill 
ridge, projected upon the (z, L) plane, follows a straight line with a slope Az/ AN given by the standard WLC model. 
When L enters the domain 20 < L < 30 the ridge bends under the barrier repulsion to become parallel to the plates. 
At the same time the hill slope becomes steeper under the combined effects of the stretching force and the confining 
plates. 

In Fig. 5 we present the results of a quantitative analysis of the terminal monomer statistics described previously in 
a qualitative way. The solid line curve gives the derivative of the chain elongation d<K z ^ '- as a function of L. When 
L stays within the interval 12 A < L < 20 A, the rate of variation per unit legnth (z(L)) stays constant as in the 
unconstrained WLC model: the elongation still behaves as an extensive quantity. We have verified that the height 
of the plateau agrees within a few tenths of a percent with the exact WLC prediction quoted above. When the 
cristallographic length of the molecular chain is going gradually from Lq = 20 A to L — L max — 30A, the elongation 
derivative with respect to L starts a rather sharp decrease towards 0. The elongation ( z(L) ) is no longer extensive 
and goes slowly to 19A ~ Lo. The dotted line of Fig. 5 gives, as a function of L, the mean square free-end fluctuations 
along the stretching force direction Az 2 (L) = ((z(L) — (z(L))) 2 ). If L < Lo, then, Az 2 (L) increases linearly with 
L and the slope stays very close to that predicted by the WLC model. When L varies from Lo to L max , Az 2 (L) 
undergoes a sharp decrease and reaches a final value ten times below the prediction of the unconstrained WLC model. 

We would like to stress that a basic physical feature of our model is the Internal Confinement ( IC), which simply 
means that all the monomers are constrained practically to stay within the space domain: < z n < Lq. We will 
find helpful, later on, to consider also models involving an External Confinement (EC): the confining potential V(r) 
is acting only upon the terminal monomer, for instance via an attached bead, so that the internal monomers are no 
longer spatially constrained. This type of confinement can be implemented by adding to the elastic density (|12|) the 
velocity dependent contribution r • V r V(r). We have calculated d<K - and A z 2 L versus L using the EC version of 
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The terminal monomer statistics. 
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FIG. 5: A quantitative analysis of the terminal-monomer statistics. The full curve and the dotted curve represent respectively 
the variations with L of the two physical quantities, first, the derivative of the elongation { z(L) ) with respect to L , di - z ^ - , 
second, the longitudinal fluctuations of the molecular-chain free end: A z 2 = ( (z(L) — (z(L)}) 2 ). One shoud note that both 
(z(L) ) and Az 2 (L) lose their extensive character when L > La- 



the WLC model with the confining potential of Fig. 2. These two quantities exhibit a decrease within the interval: 
Lo < L < L max , but they stay well above the values given in Fig. 5. For instance, when L goes from 25 A to 
L max = 30A, d{ z[ d l ] } decreases only from 0.4 to 0.3 and Az 2 (L) from 1 to 0.6. 



2. The internal monomer statistics. IC versus EC model 



All the curves we have presented so far concern the statistics properties of the terminal monomer of molecules 
having their crystallographic length within the range 10^4 < L < L max . In the case L = L max , it is clearly of interest 
to investigate the statistics of internal chain-segments with monomer numbers n within the range TV — AA" < n < N, 
where A" = L max /b and A A" = Lo/b. The internal monomer statistics properties can be obtained directly within the 
WLC model, but at the price of non trivial modifications of the Mathematica codes used to get the results displayed 
in Fig. 3, 4, 5. This would lead to developments falling outside the scope of the present paper. Nevertheless, a rather 
good physical understanding of what is going on near the confining barrier can be obtained within the " Gaussian 
"model, often used to discuss flexible polymers. It is defined by the elastic linear density: £9 aus — _L (j-) 2 -f V(r); 
a -1 is proportional to the rigidity of the harmonic potential acting between nearest- neighbour effective monomers. If 
we indentify V^(r) with the potential V(z) of Fig. 1, the internal monomer probability distribution, (x n ,y n , z n ), 
factorizes into three independent distributions relative to each component, (we remind that n is the effective nonomer 
number related to the coarse-grained s variable by s = nb.) In the continuous limit b/a <C 1 , our simple statistical 
model is easily solved by exploiting the analogy with a QM problem (see Appendix); it involves the following simple 
Hamiltonian: 

H z = H z0 + V(z) , H x0 = ~-^. (32) 
The internal-monomer probability distributions relative to the IC model read as follows: 

Pic{n,z n ) = Jdz N exp(fz N )(z N \exp(-b(N-n)H z ^\z n )x 

(z n \exp(-bnH z )\z ). (33) 
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FIG. 6: Differences between the internal monomer statistics depending upon the type of the Gaussian confining model. In the 
internal-confining (IC) model, all the monomers feel the confining potential. This contrasts with the external-confining(EC) 
model where the confining plate, at Lo = 20 a, is acting only upon the terminal monomer. The solid and dashed curves stand 
respectively for the IC and EC probabilitiy distributions P(z, n) for Lint = n b — 25 a, 27 a, 29 a, in that order from left to right; 
the thick lines are relative to the terminal monomer. The flattening of the terminal segment - 20 % of the monomers- against 
the repulsive barrier is conspicuous in the IC model while nothing similar is visible in the EC case, where internal monomers 
are allowed to wander beyond the plate. 

On the other hand, the formula giving PEc(n, z n ) is obtained by perfoming in the r.h.s. of the above equation the 
following replacements: 

fz N — > fz N - bV(z N ) , H z — > H z0 . 

The QM problems associated with both H z o and H z can be easily solved analytically by taking for V(z) a square 
well with a depth 3> ksT . The results displayed in Fig. 6 have been obtained with model parameters leading to 
basic physical quantities as close as possible to those appearing in the previous WLC analysis: L max = Nb = 30 a; 
L int = nb = 25 a, 27 a, 29a; L Q = 2/3L max = 20a; a / = 0.775. This last condition leads to : (z(L))/L = 0.775 in 
absence of confinement. 

In Fig. 6 the solid and dashed curves refer respectively to the IC and EC model. The two thick lines are relative to 
the terminal-monomer distributions. The z-probability distributions for internal monomers with L; mt — 25 a, 27 a, 29 a 
-appearing in that order from left to right on Fig. 6 - exhibit very striking differences between the IC and the EC 
model. The IC curves suggest that the terminal chain-segment, involving 20% of the internal monomers, is subject 
to a progressive flattening against the repulsive barrier. In contrast no such effect is observed with the EC curves: 
the EC internal monomers wander rather freely across the repulsive barrier while the terminal monomer is practically 
stuck to it, as it should . 

III. PERSPECTIVES AND POSSIBLE EXTENSIONS. 

We would like to suggest some possible applications or extensions of the work presented in this paper. 

• The tranverse fluctuations of a single dsDNA molecule attached at both extremities to a plane surface have been 
observed recently 26]. The experimental data have been analysed with succes by a dynamical model derived 
from the WLC model, ignoring the barrier effect of the surface. We believe it is of interest to understand this 
remarkable result within the spatially constrained WLC model, using a closely related configuration, more easy 
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to handle than the actual one. Only one extremity of dsDNA segment is assumed to be anchored to the surface 
while a pulling force parallel to the surface is acting upon the free end. The Boltzmann factor associated with 
the surface barrier effect will be described by the rounded-off step function introduced in Section II. 

The dsDNA configurations studied in Section II were not fully realistic since we have ignored the spatial ob- 
struction coming from the magnetic (or optical ) tweezers. In practice the dsDNA is attached near the bottom 
of the spherical bead. In recent micro-manipulation experiments |27t the tendancy is to use relatively short 
dsDNA segments with 2000 base pairs, corresponding to L — 680 nm. For a stretching force F — 0.3 pN, the 
molecule explores the surface around its anchoring point over a distance of the order of the root mean square 
tranverse fluctuation: yj (x 2 ) — yj (z) (kg T)/F ~ 1.55 A — 85 nm, which is about ten times smaller than 
the bead diameter. The DNA molecule does not really feel the curvature of the bead surface; therefore the 
bead can be approximated by its tangent plane. Ignoring the longitudinal fluctuations, we substitute to the 
obstructive bead a plane barrier lying at a distance I = (z) from the anchoring plate. To compute the bead 
obstruction correction, we suggest to follow a method similar to that used in Section II, but with one difference: 
instead of computing the elongation in terms of the pulling force, we shall rather get the force as a function of 
the elongation. Let Fq be the force associated with the given elongation I = (z), when one ignores the bead 
obstruction. The next step will be to compute the elongation in presence of the plate, subsitute for the bead, 
starting from a pulling force F < Fo and letting it grow until it reaches the value F\ where the elongation is 
back to its inital value I. The bead obstruction correction to the force associated with a given elongation / will 
then be estimated as A F = F\ — Fq. An improvement of the precision of the computations, with respect to that 
obtained in Section II, will be necessary. This can be achieved by taking a smaller link length b, say b — A/30 . 



In references |9 Ul0lllll | the WLC model has been generalized to a Rod Like Chain (RLC) Model, involving both 
bending and twisting rigidities. This makes possible the study of supercoiled dsDNA entropic elasticity below 
the denaturation threshlold. One can readily modify the RLC model in order to incorporate spatial constraints. 
The recurrence relation for the partition function Z n (z n ,6 n ,K) relative to a supercoiled DNA molecule, with 
a given torque T = ksTn acting upon its free end, is obtained by performing in the r.h.s. of the recurrence 
relation (|27|) the following replacement: 

1~WLc(dn+i,Qn, f) — > T~RLc(9n+i, n , f, - k 2 ) , where T RLC is given explicitly in ref. 11]. The anchoring-plate 
barrier is expected to have significant effects upon the so-called "hat curves", giving, for a fixed force, the 
relative elongation versus the supercoiling reduced parameter a. Let us take the "low" force case, F ~ 0.1 pN, 
where the RLC model "hat curve" dips steeply into the negative z region when \a\ > 0.03. This effect is 
attributed to the creation of plectonem structures which are allowed to wander in the z < half plane, because 
of the vanishing of their stretching-potential energy. Therefore, we can expect important modifications once the 
spatial constraints, which forbid the z < region, are incorporated. 



APPENDIX A: TRANSFER MATRIX VERSUS HAMILTONIAN METHODS FOR DS-DNA SUBJECT 

TO SPATIAL CONSTRAINTS. 



1. The auxiliary variable method. 

Our starting point, as in section I, is the following partition function: 



Z= I V [rjexp I - / £(s)ds I , (Al) 



£{s)=£n(r 2 )+ l -Av 2 -tv + Viv), (A2) 

with r = 4^ and r = Jjs(r)- The variable s with < s < L results from a coarse graining of the molecular chain. We 
have seen that, for a suitable choice of £o(r 2 ), the variable s coincides, to an arbitrary precision, with the arc-length 
s of a rectifiable curve. 

To compute the partition function, the first step involves the discretization of the variable s: s n = rib. The 
molecular chain is then represented by N elementary links or effective monomers with N = L/b. Assuming that 
the effective monomer length b is much smaller than the persistence length A, we replace the derivative by finite 
differences: 
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The partition function is then written as a multiple integral of a product of TV — 1 functions of the monomer-coordinates 



, JV-1 JV 



Z= f H d 3 r„{flcxp(-6^ sc (n))}Z m (r 1 ,r ). (A4) 

n=l n=2 

We are going to use a standard field-theory trick to eliminate high order derivatives whithout having to plug by hand 
Dirac delta functions into the functional integral. This involves the introduction of an auxiliary dynamical variable 
u. The starting point is the integral identity, written for fixed n: 



expl-^frO 2 }^^)- 1 



/°° b A 

d 3 u n exp{ — -((u Il ^4 + ^r Il ) 2 + r 2 )}, (A5) 
-oo ^ 



where C(A) = f^™ d 3 u„ exp{— ^4 (u„ + iA f„) 2 } = (^p-) * • Expanding the argument of the exponential we get the 
final identity: 

cxp{-^(r„) 2 } = C(A)- 1 J d 3 u„exp{-6 (^u 2 + i u„ • r B J } . (A6) 

Ignoring the overall factor C(A)~ N , we arrive in this way to an expression of £(s) involving, now, the two variables 
r and u: 

£(s) -» £ o(r) 2 +l Au-r+^jU 2 - f • r + V(r) . (A7) 
We perform the following integration by part: 

/ ds u • f = [u • r ]q - / ds (f • u) . (A8) 
Jo Jo 

Ignoring the "surface term", which can be estimated to be of the order of 1/N, we obtain the new expression of the 
clastic-energy density: 

£ new (s) = £ (r 2 ) - i u • r - f • f + V(v) + ^ u 2 . (A9) 

We have now a "kinetic term", which depends only upon first-order derivatives, so that the construction of the 
quantum Hamiltonian from the functional integral is now a relatively straightforward affair. However, a price has 
to be paid: some extra-work will be required to elucidate the physical meaning of the auxiliary dynamical variable 
u. We must now write the discrete version of £ new (s). Using the same kind of prescription as in Section I, one gets 
easily: 

£2S2(n) - £o( (r " ~p~ l)2 ) - ' K - u„_i) + V(r n ) + ^ u 2 . (A10) 

We have omitted the constant force f contribution since it can be easily added up at the end. We note that the 
"kinetic term" 

K (r n - r„_i,u„ - u„_i) = £o(— — j^— — ) - p( r n - r„_i) • (u„ - u„_i), (All) 

is symmetric under the exchanges: r n <-» r n _! , u„ <-> u„_i, while the "potential" contribution U(r n ,u n ) = 
V(r n ) + ya u n is clearly not. We can make £2i^( n ) symmetric by simply doing the replacement U(r n ,u n ) — > 
(U(r n , u„) + U(r n , u„_i)) /2. We not only remain within the leeway involved in any discretization procedure, but we 
also obtain an improvement since a trapezoidal integration is better than a rectangular one. 

The transfer matrix associated with the symmetrized energy density ^^{n + 1) reads as follows : 

T"«"(r„ +1 ,u„ +1 |r„,u„) = exp-6^?»(n+l) 

= cxp -b {K (r n+1 -r„,u„+i - u„) + U(r n+1 , u„+i)/2 + U(r n , u„)/2}. (A12) 

We define the transfer operator T new associated with the transfer matrix by writing: 

(r„ +1 ,u„ +1 |f" e >„,u„) =r ne ™(r n+1 ,u n+1 |r n ,u n ). (A13) 
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In the following, X will stand for linear operators acting upon functions f(r n ,u n ) of the dynamical variables. If 
f(r n ,u n ) is an eigenf unction, then the associated eigenvalue will be written without hat: X. 
Let us isolate the "kinetic" part of T new by writing: 

f new = exp (-± U^j TS ew exp (~U\, (A14) 

(r„ + i,u„ + i|T n ™|r„,u„) = exp (~bK (r n+1 -r n ,u n+1 - u„)) , (A15) 

where U is the operator associated with the potential U(r n ,\i n ) and -£To( r n+i — r„,u Il+ i — u„) the "kinetic" elastic- 
density term given by ea. (|All|l . We note that the "kinetic" transfer operator T$ ew is invariant upon space translations: 
r — > r+a and u — * u+b. This implies that T™ et0 is diagonal within the momentum basis defined by ( r |p r ) = exp ip r -r 
and ( u|p u ) = exp ip u ■ u. Performing the involved Fourier transforms, we get readily the matrix element of Tq ew in 
the momentum basis: 

(Pr„ + 1 ,Pu„ + 1 |f " e -|p r „, Pu J = (2^) 6 <5 3 ( Pr „ +1 -p r „)<5 3 ( Pu „ + 1 -p u „)f ne -( Pr „, Pu J, 

?<T"(Pr„,PuJ = (2^) 3 exp{-6(£ ( p 2j + iPu?i .p r J}. (A16) 

The two S- functions appearing above reflect the translation invariance of the transfer operator To. We note that p u „ 
has taken the place of the velocity v„ = r„ in £o(r^). This suggests that p u can be identified with the velocity v, as we 
are going to prove later on. We have now all what we need to write To as the exponential of an Hamiltonian operator 
Hq symmetric in Fourier space, but non-Hermitian. Returning to the ordinary space, we introduce the conjugate 
momentum operators relative to the dynamical variables r and u: 

p r = -«V r , p u = -«V U • (A17) 

Let us remind that the transformation /(p r ) — > Pr/(Pr) 7 written in the Fourier space, reads in ordinary space 
/(r) — ► p r /(r) ( similarly for the variable u); so we can write T^ ew as the differential operator exp(— b Hq cw ), where 
fjnew j g gj verLj U p t an irrelevant additive constant, by the following expression: 

£ " etu =£o(p u ) + *Pu-p r . (A18) 
The transfer operator T new , defined in eq. l|A14|) . reads then as follows: 

f new = exp f-~ U^j exp(-bH™ w ) exp L± U^j . (A19) 

Using the Cambell-HausDorff formula one can derive the following expression for T new [2j|: 

f new = exp -b {hS"" + U + ^(2 [ [U, H 7 iew ],U}} - [ [U, H$ ew ],H£ ew ]) + 0(b 3 )^j . (A20) 

The total Hamiltonian H new is readily obtained as the lowest order term with respect to b: 

H new = H^ w + U = £ (pD + i p u ■ Pr + V(r) + ^ u 2 . (A21) 

Another interest of the above formula is to allow a comparison of the symmetric transfert operator T new connecting 
adjacent effective monomers along the discretized chain with the exact evolution operator exp —bH new associated with 
a jump As — b along the original continuous chain. The formula 1A20|1 shows clearly that the absolute uncertanity 
introduced by using T new instead of the exact evolution operator is of the order of 6 3 . The use of an unsymmetrized 
transfer matrix would have led to an absolute uncertainty of the order b 2 . 



2. Physical interpretation of the auxiliary variable u : the final form of the Hamiltonian. 

In the above derivation of H new , the auxiliary variable u and its conjugate momentum operator p u have been intro- 
duced in a rather formal way. We would like to show that the operator p u has a well defined physical interpretation. 
The fact that it replaces r in going from £ new to H, suggests that these two quantities are identical. 
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To get an explicit proof, it is convenient to introduce the "Euclidian Heisenberg" operator relative to the monomer 
coordinate along the chain: 

r(s) = exp(s H) r exp(-s H). (A22) 

These operators enter in the computation of the correlation function of the two-momomer coordinates which is given 
by the standard formula [23l |: 

r Tr((exp-LH)f i (s 1 )f j (s 2 fj /Tr (exp -L H\ . (A23) 

The two-monomer correlation function relative to the velocities is obtained by taking the partial derivatives of 
C (ri(si),rj(s2)) with respect to si and s 2 - We have then to compute the s-derivative of f(s): 

= dr(s) = n exp (_ s h) = pjs) . (A24) 

as 

Going to the limit s — > the operator p u is just the Euclidian velocity operator v = [H, r]. It is clear that the above 
computation can be applied to more general situations where fj(s 2 ) is replaced by any physical operator. Let us 
choose, for instance, the Unit operator 1. One obtains readily a simple physical confirmation of the p u <^ v identity 
through the following relations between thermal averages: 

(f (*)> = <v( S )> = <p u ( S )>. (A25) 

In all the explicit computations to be performed later on , we shall use a functional basis diagonal with respect to the 
velocity vector operator v$ = v where the real vector v stands for the three corresponding eigenvalues. We can use 
a more familiar language by working with the Fourier transforms of the wave functions: ^(u) = J d 3 v exp(iu • v) $(v). 
The two linear transformations ^(u) — * p u vp(u) and ^(u) — > u'I'(u) read respectively in Fourier space: $(v) — > 
v$(v) and $(v) — » i V v $(v). We get the final expression for the "physical" Hamiltonian H by performing upon 
jjnew re pi acemen t p u _ > v and u->i V v : 

H = -^V v 2 + £ (v 2 ) + v (V r - f) + U(r) . (A26) 

As noted before, H is not an Hermitian operator, but if V(r) is an even function of r, it is self-adjoint with respect 
to the functional scalar product: (<!>! | <J> 2 ) = / rf 3 v rf 3 r $i (v, — r) $ 2 (v, r). 



3. The Spatially Constrained WLC Hamiltonian. 



In order to apply the present formalism to the WLC model with spatial constraints, we have to take an adequate 
functional form for £o(v 2 ). As it was shown previously in section I, the solution is rather simple: one introduces the 
small length 8b such that 8b/b 1 and onchooses for £q(v 2 ) the following expression: 

^) = b^L. (A27) 

Introducing the above choice of £o(v 2 ) in the final expression for the Hamiltonian H given by equation (|A26Jl . we get 
the Hamiltonian H(Sb) : 

S{Sb) = " Ja Wv2 + b ^W~ + v ' (Vr " f) + v{r) ■ (A28) 

Performing the trace upon v = |v| upon the transfer operator exp— bH(v, 5b) 1 only the value v = 1 does contribute 
in the limit 5 b — ► 0, so that v coincides with the unitary tangent vector t appearing in the WLC model. We arrive in 
this way to the WLC Hamiltonian adequate for spatially constrained dsDNA [TlllllllllllllllllllllllMElIll: 

Hscwlc = -7^V t 2 - f • t + t • V r + V(r) = H W lc({) + t ■ V r + V(r) . (A29) 

If V(r) = 0, it is easily seen that Hscwlc is diagonal in the momentum space. Perfoming the trace with respect 
to r upon the exact transfer matrix is equivalent, within the momentum basis, to take the limit p r — > 0. In this way, 
the standard WLC model Hamiltonian is immediately recovered pcj . 
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4. The transfer matrix deduced from the SCWLC Hamiltonian. 

It is of interest to compare the symmetric transfer matrix obtained from the Hamitonian H$cw lc , by an equation 
similar to (A19), to the one derived in Section I, directly from the the partition functional integral of ea. (|15fl . 
To proceed, we write the Hamiltonian H$cwlc of ea. (|A29|l under the following compact form: 

Hscwlc = H WLC {f -iPr) + V(r) . (A30) 

Using ea.(' IA20|) . the exact transfer operator exp(— bH SCWLC ) can be written, up to corrections of the order of 6 3 , as 
follows : 

f scwlc = exp ^-fry^ CX P (~bHwLc(f - «Pr)) exp^-6^ ■ ( A31 ) 

We can write the recurrence relation obeyed by the partition function Z n (r n , t„) using the transfer matrix associated 
with T scwlc'- 

Wn+i.WO - J d 3 r n J d? 0(t») exp (-6 Y^±ll + V ^ ) 

I TwLc(f — i Pr)| r n> u n) Z n (r n ,t n ) . (A32) 

Let us, first, evaluate the matrix clement of Twlc(£ — *Pr)> which is diagonal in the momentum basis. The corre- 
sponding diagonal element is just the WLC-model transfer matrix given in equation l|25[) for the case of a complex 
force f — ip r . We get in this way an intermediate expression of the TyyLc(f ~ ip r ) matrix element, within the 
coordinate space, written as a Fourier integral: 

(r n+ i,t n+ i \T WLC {f - «p r )|r n ,t„) = J ^ ^*" 3 exp i p r ■ (r n - r n+ i + b tn+1 ^ r tn ^ x 

(t n+l |f WiC (f)|t n ). (A33) 
The integration over p r gives a Dirac (5-function which is similar to that introduced by hand in Section I: 

(r„ +1 ,t„ +1 | f WLC {f - ipr)|r„,t„) = 5 3 (r n+1 - r„ - b tn + l + tn ^j ( t?1+1 | T WLC {f)\t n ) . (A34) 

Performing the integration over r„ leads immediately to the partition-function recurrence relation: 

Z n+1 (r n+1 ,t n+1 ) = exp V(r n+ i)j J d 2 Q(t n ) exp{-^ V (r n+1 - ^{t n +i + 1„)^ } 

(A35) 

5. Relation berween the symmetric and asymmetric transfer matrices. 

In principle, one can build a numerical iteration procedure from the above recurrence relation. It will turn out to 
be more complex and physically less transparent than the one which is actually used in the present paper. We are 
going to describe the simplifications to be made on ea. (|A35(l in order to recover the recurrence relation Ij26|l derived 
directly in Section I, by starting from the Boltzmann functional integral. 

a. Basic assumptions. In order to justify this procedure, we shall assume that the rates of variation with |r|, 
both of the potential V(r) and the starting partition function Zo(r,t), are at most of the order of the inverse of 
the persistence length A, which is assumed to be at least ten times larger than the length b of the elementary link. 
More precisely we shall impose the two constraints: Zq(t, t) _1 V r Zo(r, t) ~ A^ 1 and y(r) -1 W r V(r) ~ A^ 1 . In 
the computations presented in Section II, these constraints were indeed satisfied. One can prove by recurrence that 
these estimates hold also for Z n (r n , t n ).This last point can be verified by looking at the the probabilitity distributions 
P(z, L) plotted on Fig. 5, where A is taken as the length unit. This will allow us to make two important simplifications: 
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b. Simplification 1. The two Boltzmann factors involving the potentiel V(r) can be replaced by a single one, 
namely: exp — b U(r„ +1 ). For the confining potential V(z) used in this paper, we have found that the corresponding 
Boltzmann factor corrections are about 2%; they are localized near the two plate barriers in narrow regions of width 
~ A. Their overall effect is expected to be somewhat smaller, since the width of the potential well is > 20A. 

c. Simplification 2. The second simplification is a change of the monomer coordinate in the partition function 
Z n appearing in the r.h.s of eq. (|A35(I : 

v h n+1 — r n+ i — |(t n +i + t„) — > r° +1 = r„_|_i — bt n+ i. It is to convenient rewrite r^ +1 and r b n+l in terms of 
At„ = (t„ + i — t„): 

r° +1 =r n+1 -6t n -6At„ and r£ +1 = r n+1 - bt n - b/2 At n . (A36) 

We get expressions for the probability distributions P a (r n+ i) (in the limit f — > 0) by integration over t„ + i of the two 
sides of eq.J2Sl and then exchange of the integration order: 

P Q (r, i+1 ) = exp(-bV{r n+1 )) J d 2 fi(t„) J d 2 fi(t n+1 ) x 

exp (-^^n) z n{r n +i -bt n - bAt n , t„). 

We now make a first order Taylor expansion of Z n with respect to 6At„ and perform the Gaussian average over At„ 
using t n+ i as integration variable. We arrive in this way to a formula useful for our purpose: 



P a (r n+1 ) =exp(~bV(r n+1 )) J d 2 fi(t n ) ^1 - h — V r „ ■ t r ^j Z n {v n+1 - bt n , t„). 



(A37) 



The formula for P b (r n+ i) is obtained by mutipliying by 1/2 the b 2 /A term, as it follows immediately from ea. l|A36l ) 
Using the above basic assumptions, one finds that the difference P°(r n+ i) — P b (r n+ i) is ~ | (b/A) 2 and this reflects 
the lack of symmetry of the transfer matrix used in Section I. By implementing the simplifications 1. and 2. in 
equation l|A35(l . we recover immediately the partition function recurrence relation derived in Section I fea. (|26|l : 

Z n+1 (r n+1 ,t n+1 ) = exp(-6U(r„+i)) / d 3 t n (t n+1 \f WL c({)\tn)Z n (r n+1 - 6t n+ i,t„) . (A38) 
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